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ABSTRACT 

From WASP photometry and SOPHIE radial velocities we report the discovery of WASP-40b 
(HAT-P-27b), a 0.6 Mj up planet that transits its 12 th magnitude host star every 3.04 days. The 
host star is of late G-type or early K-type and likely has a metallicity greater than solar ([Fe/H] 
= 0.14 ± 0.11). The planet's mass and radius are typical of the known hot Jupiters, thus adding 
another system to the apparent pileup of transiting planets with periods near 3-4 days. Our 
parameters match those of the recent HATnet announcement of the same planet, thus giving 
confidence in the techniques used. We report a possible indication of stellar activity in the host 
star. 
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Introduction 

While the Kepler mission is currently produc- 
ing the most c andidates for trans iting extraso- 
Borucki et al.ll2010l ) 



lar planets (e.g 
based transit-search 



the ground- 
continue to find 



programs 

more planets around stars at brighter magnitudes 
than those found in the space missions. Of these, 
Hun garian Automated Telescope Network (HAT 



net: 



Bakos et al.ll2004l) and Wide Angl e Search for 



Planets fWASP: iPollacco et all 120061) have been 



the most successful. Both projects are based on ar- 
rays of 200 mm f/1.8 lenses backed by CCDs, with 
the biggest difference being that HATnet operates 
at several longitudes while WASP consists of one 
station in each hemisphere. The two projects look 
at overlapping regions of sky, which has led to 
some near-simultan eous discoveries, such as the 
planet WAS P- lib (I West et all l2009h also being 
HAT-P-lOb (jBakos et al.ll2009t ). Reporting of such 
independent discoveries gives important informa- 
tion on the reliability of the respective techniques 
and on the completeness of the transit surveys. 
Recently HATnet announced the planet HAT- 
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P-27b (jBekv et all l201ll) . a hot Jupiter in a 
3 day orbit around a my — 12.2 star. This 
planet had been independently discovered by the 
WAS P project and assi gned the name WASP- 



40b (jHellier et al.l 120111 ). We report here on 
the discovery of WASP-40b made using data 
from SuperWASP-North and WASP-South com- 
bined, together with radial velocities from the SO- 
PHIE spectrograph at the Observatoire de Haute- 
Provence (OHP) observatory. 

2. Observations 

We observed WASP-40, an -KO-type star lo- 
cated in Virgo, with the SuperWASP-North and 
WASP-South cameras durin g the three seasons of 
2008- 2010. A transit search ([Collier Cameron et al 
2006) of the resulting 30,260 photometric measure- 
ments found a strong 3.04 day periodicity. The 
discovery light curve is displayed in Figure [1^, 
folded on this period. 

Using the SOPHIE spec trograph mounted on 
the 1.93 m OHP tel escope (|Perruchot et al.ll2008t 
Bouchv et all l2009h . we obtained eight spectra 
of WASP-40 during 2010 April and May. The 
high-efficiency mode and slow readout were used, 
corresponding to a spectral resolving power of 
40 000 and lowest readout noise. We acquired a 
spectrum in both entrance fibers of the spectro- 
graph to allow monitoring of and correcting for 
the sky background. However, over the sequence, 
this background was low enough that we did not 
have to apply such a correction. Signal-to-noise 
ratio values range from 22 to 35, for exposure 
times of 30 to 43 minutes. Radial-velocity (RV) 
measureme nts were computed by weighted cross - 
correlation ( Baranne et al. 1996 ; Pepe et aLll2005l ) 
with a numerical G2 spectral template. To ac- 
count for systematic effe cts associated with the 
high-efficiency mode (e.g., Bouchv et a,l.ll2009h . we 
added an uncertainty of 10 ms^ 1 in quadrature to 
the formal errors. RV variations were detected 
with the same period found from the transits and 
with a semiamplitude of 91 ms -1 , consistent with 
a planetary mass companion. The RV measure- 
ments are listed in Table Q] and are plotted in Fig- 
ure [2] 

To test the hypothesis that the RV varia- 
tions are due to spectral line distortions caused 
by a blended eclipsing binary or starspots, we 
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Fig. 1 — Photometry of WASP-40, with the best- 
fitting transit model superimposed, (a) WASP 
discovery light curve, folded on the ephemeris of 
Table HI (b) WASP data around the transit, 
binned in time with a bin width of ^11 minutes), 
(c) High-precision transit light curves from Liv- 
erpool Telescope RISE (top, red) and Euler/C2 
(bottom, blue). 
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Table 1: SOPHIE Radial- Velocity Measurements 



BJD - 2,400,000 a 


RV 




BS 




(km s" 1 ) 


(km s" 1 ) 


(km s- 1 ) 


55,301.4903 


-15.796 


0.012 


-0.035 


55,303.4817 


-15.702 


0.013 


-0.043 


55,305.4918 


-15.837 


0.012 


-0.044 


55,323.4792 


-15.858 


0.014 


-0.059 


55,324.5650 


-15.743 


0.014 


-0.025 


55,334.4374 


-15.630 


0.014 


-0.037 


55,335.5319 


-15.809 


0.018 


-0.054 


55,336.5056 


-15.717 


0.013 


-0.014 



BJD = Barycentric Julian Date (UTC). 
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Fig. 2.— (a) SOPHIE RVs with the best-fitting 
circular Keplerian orbit superimposed. The RV 
represented by an open, red circle was excluded 
from the analysis, as it was taken during transit, 
(b) Residuals of the RVs about the fit. 



performed a line-bisector analysis (|Queloz et al 
2001I ) of the SOPHIE cross-correlation functions 
The lack of correlation between bisector span and 
RV (Fig. [3J supports our conclusion that the peri- 
odic dimming of WASP-40 and its RV variations 
are due to a planet. 

To refine the system parameters we obtained 
high-precision transit photometry On 2010 June 
26 we obtained a full transit of WASP-40 with the 
RISE (rapid imager to search for exoplanets) high- 
speed CCD came ra mounted on the 2.0 m Liv- 
erpoo l Telescope ()Steele et al. 20081 iGibson et al.l 
20081) . RISE has a wideband filter of ~ 500-700 
nm, which corresponds approximately to V + R. 
We obtained 285 exposures in the 2x2 binning 
mode with an exposure time of 35 s and effectively 
no dead time. To minimize the impact of flat- 
fielding errors and to increase the duty cycle we 
defocused the telescope by 1 mm. On 2010 June 
30 we observed a full transit of WASP-40 through 
a Gunn r filter with the C2 camera on the 1.2 m 
Euler Swiss telescope, which was defocused by 0.1 
mm. The seeing ranged from 0.5 to 1.1", and air 
mass ranged from 1.2 to 1.5 over the course of the 
105 exposures. We performed differential photom- 
etry relative to several stab le reference stars, u sing 
the ULTRACAM pipeline (jPhillon et al.ll2007l) for 
the RISE data and IRAF for the Euler data. The 
resulting light curves are displayed in Figure [TJ;. 

3. Stellar parameters from spectra 

The individual SOPHIE spectra of WASP-40 
were co-added to produce a single spectrum with 
an average S/N of around 60:1. The analysis was 



perfo rmed using the methods given in lGillon et al 
(120091) . The Ha line was used to determine the ef- 
fective temperature (T e g), and the Na I D and Mg I 
b lines were used as surface gravity (log g) diagnos- 
tics. The parameters obtained from the analysis 
are listed in Table [5J The elemental abundances 
were determined from equivalent-width measure- 
ments of several clean and unblended lines. A 
value for microturbulence ( £t) was deter mined 
from Fe I using the method of lMagain (1984). The 
quoted error estimates include those given by the 
uncertainties in T c g, logg and £t, as well as the 
scatter due to measurement and atomic data un- 
certainties. 

We estimated the sky-projected stellar rotation 
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Fig. 3. — Bisector span variation with respect 
to radial velocity, with approximately the same 
scale used for each axis. The systemic velocity 
was subtracted from the radial- velocity values. We 
adopted uncertainties on the bisector spans that 
were twice the size of those on the radial velocities. 



Table 2: Stellar Parameters from Spectra 



Parameter 


Value 




5200 ± 150 K 




4.5 ± 0.2 (cgs) 


6_ 


0.9 ± 0.2 kms- 1 


v sin/ 


2.5 ± 0.9 kms" 1 


[Fe/H] 


0.14 ± 0.11 


log A(Li) 


< 0.5 


M* 


0.91 ± 0.08 M Q 


/ 2 


0.88 ± 0.22 R e 


R.A. (J2000) 


14 h 51 m 04.19 s 


Decl. (J2000) 


+05°56'50.5" 


my 


12.2 ± 0.1 


mj a 


10.63 ±0.03 




10.25 ±0.02 




10.11 ±0.02 


USNO-B1.0 


0959-0237786 


2MASS a 


14510418±0556505 



Skr utskie et a l. (2006) 



velocity (v sin/) by fitting the profiles of several 
unblended Fe I lines. For this, we used an in- 
strumental FWHM of 0.15 ± 0.01 A, determined 
from the telluric lines around 6300A and assumed 
a value for macroturbulen ce (i> ma c) of 1.0 ± 0.3 



kms" 1 (|Bruntt et al.ll2010h 



As a check of our analysis we also estimated 
the metallicity and v sin / directly from the cross- 
correlation function of the averaged S OPHIE spec 



tra, u sing the methods described in iBoisse et al 
(|2010h . These give [Fe/H] = 0.06 ±0.09 and v sin/ 
= 3.3 ± 1.0 kms , which are in agreement with 
our preceding values. 

We input our value s of T e g, logq and [Fe/H] 
into the calibrations of Torres et al. I (l2010h to ob- 
tain estimates of the stellar mass and radius (Ta- 
ble E). 

4. System parameters from RV and tran- 
sit data 

We determined the system parameters from a 
simultaneous fit to the data described in § [2] The 
transit light curve was model ed using the formula- 
tion of Mandel fe A gol (2002 ) with the assumption 
that /?p<C/?*. Limb-darkening was accounted for 
using a four-coefficient nonlinear limb-darkening 
model, using fixed coefficients (Tabled appropri- 
ate to the passbands and interpolated in effective 
temperature, surf ace gravity an d metallicity from 
the tabulations of Claret ( 2000l ). 

The simultaneous fit was performed using the 
current version of the Markov-chain Monte Carlo 



(MCMC) cod e described bylCollier Cameron et al 



(|2007l ) and iPollacco et al.l (|2008h . The transit 



light curve is parameterized by the epoch of mid- 
transit To, the orbital period P, the planet-to-star 
area ratio (Rp/R*) 2 , the approximate duration 
of the transit from initial to final contact T14, 
and the impact parameter b = a cos i/R* (the 
distance, in fractional stellar radii, of the transit 
chord from the star's center). The radial- velocity 
orbit is parameterized by the stellar reflex veloc- 
ity semiamplitude K*, the systemic velocity 7, 



Table 3: Limb-Darkening Coefficients 



Light curve (band) 


ai 


ai 


03 


04 


WASP/Euler (RC) 


0.714 


-0.651 


1.335 


-0.602 


RISE (VJ) 


0.657 


-0.626 


1.479 


-0.658 
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and -y/ecosu; and -yesinw (j Anderson et al.ll201lf) . 
where e is orbital eccentricity and u is the argu- 
ment of periastron. 

The linear scale of the system depends on 
the orbital separation a which, through Kepler's 
third law, depends on the stellar mass M*. At 
each step in the Markov chain, the latest val- 
ues of stellar density p*, effective temperature 
T e flf and metallicity [Fe/H I are input t o the em- 
pirical mass calibration of lEnoch et al. ( 2010f ) to 
ob tain M*. The shapes of the t ransit light curve 



(jSeager fc Mallen-Ornelad 120031 ) and the radial 



velocity curve constrain /?*, which combines with 
Af* to give i?*. T c ff and [Fe/H] are proposal 
parameters constrained by Gaussian priors with 
mean values and variances derived directly from 
the stellar spectra (Table [5]). 

As the planet-star area ratio is constrained by 
the measured transit depth, Rp follows from i?*. 
The planet mass Mp is calculated from the mea- 
sured value of Ki and M*; the planetary density 
p P and surface gravity loggp then follow. We also 
calculate the blackbody equilibrium temperature 
Tp ; A=0: where A is albedo, assuming efficient re- 
distribution of heat from the planet's presumed 
permanent day side to its night side. 

At each step in the MCMC procedure, model 
transit light curves and radial- velocity curves are 
computed from the proposal parameter values, 
which are perturbed from the previous values by 
a small, random amount. The x 2 statistic is used 
to judge the goodness of fit of these models to the 
data and a step is accepted if x 2 is lower than for 
the previous step. A step with higher x 2 is ac- 
cepted with a probability exp(— Ax 2 /2). To give 
proper weighting to each transit and RV data set, 
the uncertainties are scaled at the start of the 
MCMC so as to obtain a reduced x 2 °f unity. 

From an initial MCMC fit for an eccentric or- 
bit, we found e = 0.13^oio- The improvement 
in the fit resulting from the addition of -y/ecoscj 
and Y^sinw as fitting parameters is too small to 
justify ado ption of an eccentric orb it. The F-test 
approach of lLucv fc Sweeney (1971) indicates that 
there is an 84 % probability that the improvement 
in the fit could have arisen by chance if the under- 
lying orbit were circular. In the absence of con- 
clusive evidence to the contrary, we adopted the 
circular orbit model. 



One spectrum of the eight SOPHIE spectra 
was taken during the start of transit. As we did 
not fit for the Ro ssiter-McLaughlin effect (e.g., 
Queloz et al1l2000h we excluded the resulting RV 
measurement (BJD = 2,455,301.4903) from our 
analysis. 

The median values and 1 a uncertainties of the 
system parameters derived from the MCMC model 
fit are presented in Table 2J The corresponding 
transit and orbit models are superimposed on the 
transit photometry and radial velocities in Fig- 
ures [T] and [21 

For a transit to be grazing, the grazing criterion 
(|Smallev et al.ll2011l ) must be satisfied: 



X = b + R P /R i , > 1. 



For WASP-40b X = 0.9895; 



0.014- 



(1) 

Figure [4] shows 
the MCMC posterior distribution of Rp / i?* and b. 
A total of 26.4% of these data points satisfy the 
grazing criterion. U sing the odds ratio test (e.g., 
Kipping et al. 2010l ). we find a 40.5% probability 
that the system is grazing. 

5. System age 

Assuming aligned stellar-spin and planetary- 
orbit axes, the measured vsinl of WASP-40 and 
its derived stellar radius (Table [2]) indicate a rota- 
tional period of P ro t = 17.8±7.8 days. Co 
this with the B — V color of a K0 star from Grav 



(2008j), we used the relationship of iBarnesI (|2007 
to estimate a gyrochronological age of l-2+ J Gyr. 
Considering that the stellar-spin axis may not be 
in the sky plane, these are upper limits on the stel- 
lar rotation period and the gyrochronological age. 
We found no evidence for rotational modulation 
in the WASP light curves. 

We interpolated the s tellar evolutio n track s of 
Marigo et all (|2008l) and iBertelli et all (|2008l) us- 
ing from the MCMC analysis and using T c s 
and [Fe/H] from the spectral analysis (Fig. [5]). 
This suggests an age of 6 ± 5 Gyr and a mass of 
0.83 ± 0.07 M Q for WASP-40. 

6. Discussion 

As we find more exoplanets that transit their 
host stars we begin to see patterns in their dis- 
tribution. It is thus important to add new sys- 
tems to probe such patterns and to understand 
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Table 4: System Parameters from RV and Transit 
Data 



Parameter 



Value 



P 

Ti a 
c 

T14 

12 ~ ^34 
AF - i?2 /R 2 

6 

i 

Ki 

1 
e 

M* 

i?* 

log 9 

p* 

Tcff 

[Fe/H] 

M P 

i?p 

log^p 

Pp 

a 

?P,A=0 



3.0395589 ± 0.0000090 days 
HJD 2, 455, 362.31489 ± 0.00023 
0.0696 ±0.0011 days 
0.0282t™ d ° 3 fi 7 nod days 

n ni r;o+0.00041 
U.UlOZ_ QQ057 

0-866j£8iS 

85.oi±8:fg ° 

91 ± 13 ms" 1 
-15748.7 ± 1.0 ms^ 1 
(adopted) 

0.921 ± 0.034 M Q 
0.864 ±0.031 R Q 
4.529 ±0.027 (cgs) 
1.43 ±0.13 p© 
5246 ± 153 K 
0.13 ±0.11 



0.617 ±0. 



M 



.Tup 



1.038± ; 050 R Jup 
3.112 ±0.080 (cgs) 
0.54 ±0.12 p Jup 
0.03995 ± 0.00050 AU 
1177 ±42 K 



a HJD = Heliocentric Julian Date (UTC). 
b 1 a upper limit undefined as system is near- 
grazing. 




Fig. 4 — The MCMC posterior distributions of b 
and Rp/R*. The solid line indicates the position 
of the stellar limb, i.e., b + Rp/R*= 1. A total 
of 26.4% of the points lie above the line and are 
grazing solutions. 



the role of selection effects in transit surveys. 
For example, given the apparent pileup of tran- 
siting planets with periods near 3-4 days, it is 
important to investigate the shap e of the cut- 
off t oward shorter periods (e.g., Szabo fc Kiss! 
20111 ). Further, there are increasing suggestions 



of correlations between planetary radii, the ir- 
radiation of the planet and the metallicity of 
the host star (lEnoch et aL 2011 ; Anderson fc Irol 
2011. in preparation! ) . For brighter systems, the 
measurement of the Rossiter-McLaughlin effect 
allows us to build up st atistics of the align ment 
of planetary orbits (e.g., Triaud et al. 2010l ). and 
to thus test suggestions of a correlation between 
alignment and the spectral type of the host star 
dWinn et alJl20ld) . 

Since all such comparisons depend on the re- 
liability of measured parameters of the transit- 
ing systems, it is worth noting that the param- 
eters of HAT-P-27b/WASP-40b reported here are 
in good a ccord with t hose f rom the independent 
study by iBekv et"al1 ( 201 lh . Most parameters 
agree to within 1 a errors, while the impact param- 
eter, which is strongly dependent on the modeling 
of the transit light curve and the limb-darkening, 
agrees to within 2 a. This adds confidence to the 
methods used by the two projects. 

One difference in our analyses is that our radial 
velocities show more scatter about the model than 
would be expected from their formal uncertainties 
(Fig. IH]). This could be an indication of stellar ac- 
tivity. We calculated an activity index suggestive 
of an active star, logi?^ K = —4.63, directly from 
th e SOPHIE spectra using the methods described 
m iBoisse et all (|2010l ). The RV dispersion arising 
from the stellar activity of a K-type star is esti- 



mated to be ~10 ms 1 ( Santos et al. I l2000h . 



We note that the radial- velocity data in Bek v et al 
(|201l[ ). while resulting in a very similar model, 
do not show a similar scatter, which could indi- 
cate that the level of stellar activity fluctuates. 
Alternatively, if the scatter was indicative of an 
additional planet, the sig nal may not be pr esent 
in the radial velocities of IBekv et al.l (|201l[ ) due 
to limited sampling. Both the SOPHIE and Keck 
data sets are sparse, with eight and nine radial 
velocities, respectively, so more observations are 
needed to reach a conclusion. 

The research leading to these results has re- 
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tionary mass tracks ( Z = 0.026 » [ F e/H] = 0.14; 
Y = 0.30) are from iBertelli et all ([20081 ) . The 
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ages 1, 4, 7, and 10 Gyr are from iMarigo et al 
(|2008h . 
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Fig. 6. — Residual radial velocities about the best- 
fitting circular Keplerian orbit. 
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